Charge-Resonance Effect on Harmonic Generation by Symmetric Diatomic Molecular 

Ions in Intense Laser Fields 
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In present paper we develop an analytic theory for the harmonic generation of symmetric diatomic 
molecular ions beyond two-level model, emphasizing the influence of charge-resonance (CR) states 
those are strongly coupled to electromagnetic fields for large internuclear distance. With taking 
into account the continuum states that is ignored in the two-level model and become important 
for intense laser case, our model is capable to produce spectrum for the whole range of harmonic 
orders consisting of a molecular plateau due to the CR transition and an atomic-like plateau for a 
long-wavelength excitation, and in good agreement with numerical results from directly solution of 
the Schrodinger equation. Our theory also identifies the crucial role of the CR states in the fine 
structure of harmonic spectrum and shows that the harmonic generation in molecular system can 
be effectively controlled by adjusting the internuclear distance. 

PACS numbers: 42.65.Ky, 32.80.Rm 
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I. INTRODUCTION 

In recent years, high order harmonic generation 
(HOHG) from atom and molecule has been the subject 
of numerous experimental and theoretical studies, mainly 
because the emission of high-order harmonics is a promis- 
ing method to produce coherent x rays and attosecond 
pulses, and an effective ways to study the internal con- 
struction of atom or molecule^ Ej Ell li l [HIE H HH ITol 

IH 13 13 El El El III 13 Ellallll^ 

Compared to atom case, the HOHG spectrum of 
molecule systems demonstrates novel properties, e.g., 
peak splitting and sideband peaks, Rabi or Mollow 
triplets effect, and combination of the atomic-like and 
molecular-like plateau, molecular alignment effect, and 
two-center interference, to name only a fewpij. 124 . 12M 
|2(| . Symmetric molecular system such as have pairs 
of electronic states known as the charge-resonance (CR) 
states, (i.e.,l<7 s ,lc7 M .) It is known that the CR states 
are strongly coupled to electromagnetic fields at large in- 
ternuclear separation R. Bandrauk and co-worker first 
pointed out the importance of these CR states as sources 
of highly nonlinear laser- induced effects in molecules |23|. 
Ivanov, Corkum examined the possibility of using the 
CR states to produce and coherently control HOHG[2^|. 
Zuo, Chelkowski, and Bandrauk have showed that sym- 
metric molecular ions in general produce more efficient 
harmonic generation than atoms, especially at large R, 
due to these CR states [25|,|26|. However, in the all above 
discussions, the theoretical analysis is based on a two- 
level model where the continuum (ionization) states are 
ignored completely. Since the ground-continuum cou- 
pling and ionization become important for the strong 
laser fields of 10 14 W/cm 2 and above, a fully understand- 
ing of the whole structure of HOHG and the effect of CR 
states in molecule systems requires to extend to consider 
the continuum states and ionization process. 



In present paper we develop an analytic theory on the 
harmonic generation of symmetric diatomic molecular 
ions beyond the two-level model. Our model is capa- 
ble to produce harmonic structure for the whole range 
of harmonic order, including the molecular plateau due 
to CR transition and the atomic-like plateau for a long- 
wavelength excitation, and agrees well with the numer- 
ical results from directly solving the Schrodinger equa- 
tion. Our theory identifies the role of the CR states in 
the fine structure of harmonic spectrum and shows that 
the harmonic generation in molecular system can be ef- 
fectively controlled through CR states by adjusting the 
internuclear distance. 

Our paper is organized as follows. In Sec. II we present 
our analytical theory. Our analysis is divided into two 
cases: near-resonance region of intermediate internuclear 
distance and the strong-coupling region of large internu- 
clear distance. We will derive the time-dependent am- 
plitudes of the ground, excited and continuum states, 
then calculate dipolar moments and their Fourier trans- 
formation. The analytic expressions of the amplitudes 
of HOHG for the whole range of harmonic order is 
given in this section. Sec. Ill presents our numerical re- 
sults. Our theory is applied to ID symmetric diatomic 
molecule model, making analysis on the structure of 
HOHG and compare with the results from directly solv- 
ing Schrodinger equation. Sec.V is our conclusion. 



II. ANALYTIC THEORY 

We consider a diatomic molecule ion in single-electron 
approximation under the influence of a linear polarized 
laser field A(t) = (E cos(wt), 0, 0), where Eu> is the abso- 
lute amplitude of the external electric field and lo is the 
frequency of the external field. The Hamiltonian of the 



2 



model molecule studied here is 

H(t) = (p-A(t)) 2 /2 + V(r), 

where p is the canonical momentum and V(r) is the bind- 
ing potential. 

Under strong- field conditions [53, 13 > ^ * s reasonable 
to assume that, (a) Except the ground state |0) and the 
first excited state|l), the contribution from other bound 
states can be neglected; (b) The depletion of the ground 
state and the first excited state is small; (c) In the contin- 
uum, the electron can be treated as a free particle moving 
in the electric field without considering Coulomb poten- 
tial. Then, the time-dependent wave functions can be 
expanded as 



where S(p,t,t') = J* [(p - A(t")) 2 /2 - E ] dt". D (t) 
denotes the transition back to the ground state, and 
D\{t) denotes the transition back to the first excited 
state. Each of Do(t) and D\{t) includes two differ- 
ent terms, which denote two different moments, respec- 
tively. For D (t), they are |0) -> |p) -> |0) (denoted 
by D° (t)) and |1) -f \p) -> |0> (denoted by Dj(i)), 
which denote the transition from the ground state to the 
continuum state, then back to the ground state(L>o(t)), 
or from the first excited state to the continuum state, 
then back to the ground state(Z?Q(i)); for Z>i(t),they are 



10) 



|1) (denoted by D°(t)) and |1) 



m)) 



-iEnt 



a(t)\0)+b(t)\l) 



dpcp(t)|p>], (l) 



where — Eq is the ionization potential of the ground state, 
a(t) is the ground-state amplitude, b(t) is the first excited 
state amplitude and c p (t) are the amplitudes of the corre- 
sponding continuum states. While the ionization is weak, 
by neglecting the depletion of the ground and the first 
excited state and setting |a(i)| 2 + |6(i)| 2 = 1, the formu- 
lation of a(t) and b(t) can be obtained by a two-level ap- 
proximation. We have factor out here free oscillations of 
the ground-state amplitude by the bare frequency — Eq. 
The Schrodinger equation for c p (t) reads as 

c p (t) = iA(t)p[a(t)(p\0)+b(t)(p\l)} + [(p-A(t)) 2 /2-E } Cl 

(2) 

(2) can be solved exactly and c p (t) can be written in the 
closed form 



c p (t) = i / di'A(i')p[a(t')(p|0> + 6(t')(p|l) 



(3) 



xe 



-tf t t l [(p-A(t")) 2 /2-E ]dt" 



IP) - |1> 

(denoted by D{(t)), which denote the transition from the 
ground state to the continuum state, then back to first 
excited state(D"(t)), or from the first excited state to the 
continuum state, then back to first excited state(D}(t)). 
We should discuss these moments at intermediate R and 
large R, respectively. 



A. Intermediate R with (Ei — Eq)/uj ~ 1 

For a two-level system, the time-dependent wave func- 
tion can be written as \tp(t))=a(t)\0) +b(t)\l), where |0) 
and |1) are the ground-state and the first excited state 
of the unperturbed Hamiltonian Hq = p 2 + V(r) with 
H \0) =E \0) andff |l) = Ei\l). While (E x -E )/uj~ 
1, in the rotating- wave approximation, the solution for 
Wa(t) and b(t) can be written as 

f o(t) = e -^ot- l / t [A 2 (t')/2]dt' (aiei Q t /2 +a2e ^Qt/2y^t/2^ 
1 6(t) = e -iEit-iJ£[A. a (t')/2]dt' Q, ie iQt/2 _|_ j, 2 g-<Qt/2^ e iet/2 ) 

and by the Bessel function formula[2^ ,we obtain finally 



Using Eq. (1) and (3), the mechanical momentum p — 
A(t) component of the time-dependent dipole moment is 



I &(*) = E« MEySio^he-^ + b 2 e- iBat ), 



{ <*(*) = EnZZZ Jn{E 2 /^){ ai e~^ + a 2 e 

b 2 



(7) 



T>( t ) = (Tp(t)\p-A(t)m)). 



(4) 



Neglecting the term (ip(t)\A(£)\ijj(t))=A(t), which only 
includes the fundamental frequency u>, and the contribu- 
tion from C-C part and considering only the transitions 
back to the ground state and the first excited state, we 
obtain D(i) = D (t) + D x (t), where 

D (t) = i J dp J dt'a*(t)(0\p\p)A(t')-p[a(t')(p\0) 

+6(t')(p|l>] x e -^(p. i '*') ! 



where Euj is the absolute amplitude of the external elec- 
tric field andAi = £ 2 /4 + 2nw-Q/2 + £/2, A 2 = E 2 /4 + 
2nu; + Q/2 + Z/2,B 1 = E 1 +E 2 /4 + 2mu- Q/2- f/2- E , 
B 2 = E 1 +E 2 /4: + 2nuj + Q/2-£_/2-Eo, where £=E 1 - 
Eo-cj; Q=^\Q R \ 2 + e, Qr=ED, £>=(0|p|l) is the ma- 
trix element of the electric dipole moment, and a\,a 2 , b\ 
and b 2 are constants of integration which are determined 
from the initial conditions |^(0))=a(0)|0) + 6(0) |1>: 



(5) 



D 1 (t) = i J dpj^ dt'6*(t)(l|p|p)A(* , ).p[a(0<p|Q) 

+6(t')(p|l)] x e -' iS (p^.*') 



oi - [(Q + 0«(0) + Qh6(0)]/(2Q) 
a 2 = [(Q - C)a(O) - Q R b(0)]/(2Q) 

6 1 = [03 - 06(0) + Q* R a(0)]/(2Q) 

6 2 = [(Q + e)6(0)-Q|ja(0)]/(2Q), 



(8) 



(6) 



where Q* R is the conjugated element of Qr. For the con- 
venience of latter calculation, we have multiplied a(t) and 
b{t) by a phasic factor e iEot . Substituting (8) into (5) and 
(6), then it follows: 
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1. |0> — |i>> — |0> 



The Fourier component of D®(t) reads as 



For Do(t), above all we consider the moment |0) 
|P> -» |0) , i.e., 



n,m,ni,mi-+oo 

-i I dpEir(0\p)pp x (p\0) ^ J n J m 

n,m,rii,mi--oo 



alaiS [(m + mi - 1 + 2(n — n\))u) — u)'] 



D° (t)=i J dpj^ diV(i)<0|p|p>A(t')-pa(i')<p|' 



-iS(p,t,t') 



(9) 



For e~ iS ( p '*'*'\ by the formula e -*7s*™(2wt)+*K£sin(u;t) _ 
En=-Z Jn{-kl)e- inut , we can obtain 



e -»S(p,t,t') =e -»(^ 2 /4+p 2 /2-E )(t-t') 



n,m— +oo 

E . 

— — 

-£ 2 



n,m— — oo 
^2 



7" f ^ W / P 1 ^ _^_\ -inoit -ir?i^/ 



(Ci+^'i) 








a2a 2 <5 [(to + mi — 1 + 2(n - 


- m))w 


-a/] 




(C 1 + A' 2 ) 
a^ai(5 [(to + mi + 1 + 2(n - 


- ni))w 


-a;'] 




(C 2 + ^i) 
a2ft 2 <5 [(to + To! + 1 + 2(n - 


- m))w 


-a/] 




(<?2 + A' 2 ) 
a^a 2 (5 [(to + toi — 1 + 2n — 


2ni)w - 


+ Q- 


u/] 


(C\ + A' 2 ) 
ala 2 S [(to + toi + 1 + 2n — 


2ni)u - 


+ Q- 


u/] 


(C2 + A> 2 ) 
a 2 a\8 [(to + toi — 1 + 2n — 


2ni)w 


-Q- 


c/] 


(Ci + A[) 








a 2 a\8 [(to + toi + 1 + 2n — 


2ni)w 


-Q- 





(10) 



substituting it into (9), and integrating over t', we obtain 



^ °(t)=* y dp|<oip)p( Px )<p|o) 

n,m,ni,mi=+oo , a * aie _ l(Cl +A' 1 +_D 1 -A 1 )t 



E 



n,m,ni ,mi— — oo 



-iCCi + Ai) 



+ — 77^ 77^ + " 



-i(C 2 + A[) 



ala 2 e- i ^ +A ^ +Dl - A ^ t a* 2 a 1 e~^+ A '^+ Dl - A ^ t 



-i(C 2 +A> 2 ) 



-i(Ci + A[) 



+ 



a 2 a 2 e 



-i(Ci+A' 2 ) 

-i(C 2 +A' 2 +D 1 -A 2 )t 



-i(C 2 + A' 2 ) 



-i(C 2 + A[) 



For(10),the parity of (0|p)pp x (p|0) for the integral 
over p is even, according to the property of Besscl func- 
tion, only when to + toi is even, the value of Eq. (10) for 
the integral over p is not zero. Fourier components on the 
right-hand side of (10) can be divided into two parts, odd 
Fourier components of (m + m 1 +2n — 2n\ ± l)w, and no- 
integer Fourier components of (m+mi+2n—2ni±l)uj±Q 
.It also can be seen from (10), the Rabi oscillation of the 
ground and first excited state should not contribute to 
the parity of integer harmonic, but induces the symmet- 
rical splitting of odd harmonic. The splitting separations 
around each odd harmonic all take the same value Q. If 
Q/uj is integer of odd number, then the odd harmonic 
sidebands should coincide at the even harmonic (2nu>) 
position, thus giving rise to radiation of even harmonics. 



wWp t - t (ZzE =ML\ j — j / -p* E Mi\ 

wiicie <J m — <Jm\ u , 8w )i 'Jm 1 —'Jm 1 y u , SuJ ), 

Jn=Jn(£), .7 ni =J ni (f£),and A x = E 2 /4 + 2 ni uj- 
Q/2 + f/2, A 2 = E 2 /A + 2n x uj + Q/2 + £/2, 
C 1 =mu- p 2 /2 - E 2 /A + E Q - u, C 2 = muj - p 2 /2 - 
E 2 /A + Eq + u, A\ = E 2 /A + 2nuj - Q/2 + £/2, A' 2 = 
E 2 /4 + 2nuj + Q/2 + ^/2,D 1 = p 2 /2 + E 2 /4 - E a + miw. 



2. |i) _ | P ) _ |o> 
For |1) - |p) - |0), i.e., 



Dl(t)=i J dp J dt'a*(t)(0\p\p)A(t') ■ pb(t')(p\l)e 



is(p.t,t') 
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Analogous to |0) — > |p) — ► |0), the Fourier component of The Fourier component of D®(t) reads as 
D^(t) reads as 



n,m,rii,mi— +00 

dp J B7r(l|p)pp ;E (p|0) ^ J„J m 

n,m,ni ,mi= — 00 



n,m,fii,mi- +00 

dp J B7r(0|p)pp x (p|l) J„J m 

n,m,ni,mi— — 00 



aJ&iJ [(m + mi + 2n — 2m)cj — uj'] 



+ 
+ 
+ 
+ 
+ 
+ 

+ 





(Ci + B[) 


a 2 b 2 S [(to - 


h toi + 2n — 2ni)ui — uj'] 




(Ci+B' 2 ) 


a|6i<5 [(m - 


f toi + 2 + 2n - 2n{)u) - uj'] 




(C 2 + B[) 


036215 [(to - 


f toi + 2 + 2n - 2m)cj - uj'] 




(C 2 + B 2 ) 


a{b 2 8 [(to + mi + 2n — 2rii)uj + Q — uj'] 




(Ci + B' 2 ) 


alb 2 S [(to + toi + 2 + 2n - 2ni)w + Q - uj'] 




(C 2 + B> 2 ) 


a 2 b\5 [(to + toi + 2n — 2n\)uj — Q — uj'] 




(Ci+B[) 


a 2 b 1 S [(m + m 1 + 2 + 2n - 2n{)u) - Q - uj'] 



blaid [(to + mi — 2 + 2n — 2m)w - uj'] 
(Ci + A[) 

b 2 a 2 8 [(to + toi - 2 + 2n - 2m)w - u/] 

(Ci+A' 2 ) 
6iai<5 [(to + toi + 2n — 2ni)uj — uj'] 

(C 2 + A[) 
b 2 a 2 S [(m + mi + 2n — 2n\)uj — uj'] 

(C2 + A' 2 ) 

b\a 2 5 [(to + mi - 2 + 2n - 2m)uj + Q - uj'] 
(Ci + A> 2 ) 

b\a 2 5 [(m + m x + 2n — 2ni)uu + Q — uj'] 
(C2 + A' 2 ) 

b 2 a\5 [(m + mi — 2 + 2n — 2n\)uj — Q — uj'] 
(Ci + A[) 

b 2 ai5 [(to + mi + 2n - 2n\)ui - Q - uj'] 



(C2 + A[) 



(12) 



{G 2 + B[) 



(11) 



where the definitions of C\ , C 2 , A\ , A' 2 , are the same as 
in |0> - |p) - |0). 

For(12), only when m + m-i is odd, the value of(12) for 
the integral over p is not zero. Fourier components on 
the right-hand side of (12) also can be divided into two 
where the definitions of C\ and C 2 are the same as in parts, odd Fourier components of (m + m 1 + 2n- 2n x )uj 



|0) -» |p) -» |0), and B[ = E x + E 2 /4 + 2nuj - Q/2 - 
e/2 - S , B 2 = £i + £ 2 /4 + 2nuj + Q/2 - f/2 - £ - 

For(ll), the parity of (0|p)ppa;(p|l) for the integral 
over p is odd, only when m + m-i is odd, the value of(ll) 
for the integral over p is not zero. Fourier components on 
the right-hand side of (11) can also be divided into two 
parts, odd Fourier components of (to + mi + 2n — 2ri\)uj 
and (m + mi + 2n — 2n\ + 2)uj, and no-integer Fourier 
components of (m + mi + 2n — 2n\)u) ± Q and (m + m x + 
2 + 2n-2 ni )uj ±Q. 



For |0) 



3. |0> -> |p> -> |1) 
IP) - |1>, i.e., 



and (m + mi + 2n — 2n\ — 2)oj, and no-integer Fourier 
components of (m + mi + 2n — 2ni)uj±Q and (m + mi — 
2 + 2n-2 ni )uj±Q. 



I |1) -> |p> -> |1> 



D\{t)=i J dp^ t dt'6*(t)(l|p|p)A(t')-p6(0(p|l)^ s(p ' M ' ) . 
The Fourier component of D\(t) reads as 

/n,m,ni,rni— +oo 
dp£7r(l|p)pp x <p|l) J2 J « J ™ 
n,m.ni,mi--co 

J b\bi6 [(m + mi - 1 + 2n - 2ni)w - uj'] 



D° 1 (t)=i dp dt'b*(t)(l\p\p)A(t') ■pa(t')(p\0)e- lSip > t ' t ">+ 



(Ci+B[) 

b 2 b 2 S [(to + mi - 1 + 2n - 2m)u) - uj'] 

(Ci + B' 2 ) 
blhS [(to + mi + 1 + 2n - 2n Y )uj - uj'] 

(C2+B[) 
b' 2 b 2 5 [(m + mi + 1 + 2n - 2m)w - uj'] 

(C 2 + s 2 ) 



(13) 
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blb 2 8 [(m + mi - 1 + 2?i - 2n 1 )uj + Q - uj'} 





V^i 


+ Sol 


^1^2<5 [{TU + JTll 


+ 14 


2n - 2ni)uj + Q - uj'] 




(C 2 


+ B'o) 


62&i(5 [(m + mi 


-1 + 


2n - 2m)w -Q-uj'} 




(Ci 


+ B[) 


b^biS [(m + uii 


+ 14 


2n — 2hi)lu — Q — uj'] 




(C 2 


+ B[) 



B. Large R and {E 1 - E ) /to ~ 



When [E\ — Eq)/uj ~ 0, assuming Ei—E and 
using the Bessel function formula e -Ksin{ut) _ 
^n(C) e_M *j the two-level approximation solu- 
tion of a(t) and b(t) can be written as: 



where the value of Ci, C2, P^ and P 2 is the same as in 
|1>-|P>-|0>. 

For(13), the parity of (l|p)pp a: (p|l) for the integral 
over p is even, only when to + mi is even, P x is not zero. 
Fourier components on the right-hand side of (13) include 
odd Fourier components of (m + mi + 2n — 2ni ± 1)uj, 
and no-integer Fourier components of (to + mi + 2n — 
2n 1 ±l)u>±Q. 

Using(10)-(13), the Fourier component of D(t) reads 

as 

t — +00 

P(o/) = / * [!>§(*) + £»?(*) + I>S(*) + D\ (t)] e— '* 

£= — 00 

= P V) + FJV) + P V) + P X V). (14) 

It can be seen from (10)-(13), no matter initially the 
system being in the ground state or the first excited state 
or the coherent superposition state[2l], the right-hand 
side of (14) only includes odd harmonic with the sym- 
metrical splitting of ±Q. 

It should be noted that the rotating-wave approxi- 
mation is applicable only when (Pi — Eq)/uj ~ 1 and 
the external field is weak. However, in case of stronger 
field where tunnelling ionization is prominent, our model 
at intermediate R is considered to be an approximative 
approach. The splitting separation can also be calcu- 
lated by the effective Rabi frequency or quasienergies 
of the Flock states of the system^! |H |3(f . But for 
some specific intermediate P, the molecule ions have rel- 
ative high ionization rate in the presence of relative weak 
field and if the external field frequency uj accords with 
(Pi — Eq)/uj ~ 1, our model can give an adequate de- 
scription of the process. 

When (Pi — Eq)/uj ~ 0, according to^Sj, the en- 
ergy separation of quasienergy (Floquet or dressed) 
states around each even harmonic accords with oj q = 
AEJ (2Q' r /uj) ~ 0, (where Q' R = ED', D' = (0|x|l), 
P is the absolute amplitude of the external field here), 
in other words, when R — > 00, all of the large Rabi split- 
tings of the odd harmonics should converge towards even 
harmonic[3J. However it can be showed that in this case 
the amplitudes of all of even harmonic should be zero, so 
at very large P, only odd harmonic can be produced. 
The point can be made clearer by the latter analysis of 
(16)-(19). 



a(t) = -e 



x e~ inut + v 



n— +00 



uj Suj 



n— — 00 
■m— +00 

(o) e J «( z 



\Qr\ 



m — — co 



uj ' Suj 



2Qi 



n— +00 



*(o) E J »( : 



-\Qr[ 



(15) 



n— — 00 



UJ 8oj ' 



m— +00 



x e 



»(o) e J ™( ; 



■\Qr\ 



m— — 00 



uj ' Suj 



whereu(O) = a(0) + jg^6(0), and v(Q) = a(0) - ^6(0). 
For convenience, we have multiplied a(t) and b(t) by a 
phasic factor e lEot , respectively. If initially the system 
is in the ground state |0), then u(0) = l,v(0) = 1; if 
initially the system is in the first excited state |1), then 
u(0) = t^T! u (0) — — ra^T> if initially the system is in the 
coherent superposition state ^|0) + ^|1), then u(0) = 

f ( 1 + T&)'<°) = ^( 1 -&)- 

Substituting (15) into (5) and (6) and considering only 
the transitions back to the ground state and the first ex- 
cited state, one can obtain four different transition mo- 
ments |0) -> |p) -> |0),|1) -> |p) -> |0) ,|0) - |p) - |1) 
,|1) — > |p) — * |1). Accordingly, the Fourier component of 
the four moments can be written as: 



Ptt 



n , m , n 1 ,712 — +00 



n,m,ni,n,2 — -00 



dp— (0|p) PPa; (p|0) E 

<5[(m — 1 + rt + ni — 712)^ — uj' 
A 

6[(m + 1 + n + rz-i — n2)w — cj'] 
P 

(_l)«i+«2| u (0)| 2 + (-l)" 2 u *(0)7j(0) 



(16) 



+ (-l) ni 7j*(0) U (0) + |7;(0)| : 
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. . n,m,n 1 ,n 2 =+ca 

p V) = -i I d p£^(o|p)p P;c (p|i> J2 



*Qr 

S[(m — 1 + n + ni — ri2)tjj — uj'] 
A 

5[(m + 1 + n + n\ — nzjuj — uj'] 
B 



(_l)m+« 2 | n (0)| 2 - (-l) n i u *(0)v{0) 



n,m,ni,ri2 — — oo 



(17) 



+ {-l) ni v*(0)u{0)-\v(0)\ 2 



Jn Jm Jn± 



n,m,m ,n 2 — +oo 

= -i i d P ±^(i|p)p P2; (p|o> Yl 



MQr\ 



S[(m — 1 + n + rii — n- 2 )uj — U) 1 ] 
A 

S[(m + 1 + n + ni — n%)u} — uj 1 ] 
B 



(_l)»i+»>|«(0)| J + (-l) n2 u*(0)v(0) 



n,m,ni,ri2 — — oo 



(18) 



(-i)»V(o)t*(o)-Ko)| s 



p 1 i (c') = - i / dp^<i| P >ppx(p|i) E 



n^m^ni ,712 — +oo 



n,m,ni,7i2 — — 00 



5[(m — 1 + n + ni — ^2)0; — a;'] 
I 

(5[(m + 1 + n + ni — ri2)w — w'] 



4) ni + n2 |w(0)| 2 - (-l)" 2 w*(0)i;(0) 



(19) 



(-l) n V(0)u(0) + |v(0)| s 



J nJ vn.J n\J "d'zi 



From the above analysis, one can show that no mat- 
ter initially the system being in the ground state |0) or 
first excited state |1), only odd harmonics should be in- 
cluded in (20); while initially the system being in the 

coherent superposition state of ^|0) + ^ |l)|2lj. the 
formula (20) should include all odd and even harmon- 
ics, which corresponds to our numerical calculation (see- 
ing Fig. 6). According to the generalized Bessel function 
formula lim u ^o[J2n+i(u, v)] = 29], while \Qr\ — > 0, 
if one or all of the values of ri\ and n 2 are odd, then 

JnA lJ ¥> £)■/«. C 1 ^, C) - 0, hence while \Q R \ « 1, 
in all case, odd harmonics should primarily come from 
(16) and (19), and even harmonics should primarily come 
from (17) and (18). However, because the initial phases 
of |0) and |1) all take zero in our numerical calcula- 
tion, hence (l|p)(p|0)=-(0|p)(p|l), and Qr/\Qr\ = -i, 
the contributions to even harmonics from (17) should be 
counteracted by the contributions from (18), which imply 
that the intensity of even harmonics as a whole should 
be weaker than that of odd one, to which the contri- 
butions from (16) and (19) shouldn't counteract each 
other, and should depend on the relative phase of the 
ground and first excited state. Since other contributions 
to even harmonics from (16) and (19) are correlative to 
m — n 2 being odd, then their intensity also should de- 
pend on the value of \Qr\- The result also is accordant 
with the ultimate case of the strict degeneracy of |0) and 
|1), where E x = E , D = (0|p|l)=0 and Q R =0. In the 
ultimate case, all even harmonic from (16)-(19) should 
disappear which spells that while R — > 00, the symmet- 
rical diatomic molecule ions should be equivalent to two 
unattached atoms. Furthermore, since in the ultimate 
case R — > 00 and (E\ — E )/uj — > 0, so uj q — > 0, the 
Floquet states theory predicts the appearance of all odd 
and even harmonics while the system initially being in 
the ground state|l3l |25(, and our model only predicts 
odd one in the case, it seems safe to conclude that the 
even harmonics coming from the large Rabi splittings of 
the odd harmonics should get weaker and weaker and 
disappear finally with the increasing R (seeing Fig. 2 and 
Fig.6(a)). 



whereJ„ — J n ( ^ , $ u ))Jm — Jm( P ^ 1 SuJ ),Jn 1 — 
■M^.fb), Jn 2 = Jn 2 (^,£), A= (m-l + ni )w- 

^- + E ,B={m + l + ni)w - ^ + E . 

From above formulae, it is easy to prove that if initially 
the system is in the ground or first excited state, only odd 
harmonic can be generated and if initially the system is 
in the coherent superposition state, both odd and even 
harmonics can be generated. 

Using(16)-(19), the Fourier component of D(t) reads 



dt[D° (t) + DUt)+D x (t) +D{(t)] e 



-iuj t 



= p V)+aV)+^V)+AV)- 



(20) 



III. NUMERICAL RESULTS 

In this section, we apply our theory to concrete model 
and present numerical results. For convenience, we 
adopt an ID symmetrical diatomic molecule model. 
The Hamiltonian of the model molecule studied here is 



-xEsm(u)t), 



2dx2 y/l + (x+a.5R) 2 ^l + O-0.5fl) 2 

where z is the effective charges, R is the internuclear 
separation, E is the absolute amplitude of the ex- 
ternal electric field here and uj the frequency of the 
external field. In the paper, we adopt atom unit, 
fi = e = m e = 1, the laser intensity varies from 
10 13 W/cm 2 to 10 14 VF/cm 2 (0.017a.u. - 0.055a.u.), the 
internuclear separation varies from intermediate R 
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Figure 1: Photon-emission spectrum of symmetrical diatomic 
molecular ions and the population of the ground state(solid 
curve) and the first excited state(dotted curve) at R = 5.2a.u., 
I — 6.7 x W 13 W/cm 2 . and u = 0.05642a.u., initially in the 
ground state: (a)and(b) Two-level calculation; (c)and(d) ID 
time-dependent exact calculation. 



(5.2a.u., and 6a. u.) to large R (16a. u.), the wavelength 
adopted is A=800nm, accordingly id = 0. 05642a. u. and 
the laser pulse contains 20 optical cycles. Numeri- 
cally the above schrodinger equation can be solved by 
operator-splitting method|22j|. 

We begin the calculation by assuming the electron ini- 
tially is in the ground state. But Ivanov et al show that 
one should consider at least two wave packets excited to 
the electronic surfaces in their program]^, which means 
initially the system should be in the coherent superposi- 
tion state|2l|. So we also perform the calculations from 
the coherent superposition state of the ground state and 
the first excited state with equally weighted populations. 



A. Failure of the two-level modle in producing the 
whole harmonic spectrum 

When R — 5.2a.u., the ionization potential of ground 
state is Eq = — 0.89830a. u. and first excited state is 
Ei = — 0.84085a. u.. The energy separation between 
Eq and E\ is AE = 0. 05745a. it., which is in the near- 
(one-photon)-resonance region. Fig.l shows the spec- 
trum from ID time-dependent calculation for the model 
molecule at R = 5.2a.u.. In Fig.l, each odd harmonic 
peak is accompanied symmetrically by two strong side- 
bands. The energy separations between the neighboring 
sidebands all have the same value, which can be explained 
in a dressed molecule state picture [33) . 



Figure 2: Photon-emission spectrum of symmetrical diatomic 
molecular ions at R — 16a.it., I — 10 14 W/cm 2 , and ui — 
0. 05642a. u., initially in the ground state:(a) two-level calcu- 
lation; (b) ID time-dependent exact calculation. 



Comparison between Fig.l (a) and Fig. 1(c) shows that 
while the external field becomes stronger and although 
the ground states and the first excited state depletion 
in Fig. 1(c) is only 0.03, there is very large difference be- 
tween ID time-dependent exact calculation and two-level 
calculation. So the continuum state effects on spectrum 
must be considered in this case. Otherwise, an obvious 
characteristic, different from atom, for the population of 
molecule states is the strong-coupling of the ground state 
and the first excited state. It shows that it is necessary 
and reasonable to consider the first excited state effect on 
harmonic generation to explain the spectrum from sym- 
metrical diatomic molecular ions. 

While R = 16a. u., the ionization potentials of the 
ground state and the first excited state are Eq — 
-0.732574a.u. and E x = -0.7325723a.u., respectively. 



The energy separation is AE — 1.7 x 10 



which is in 



the strong-coupling region. When / = 10 14 W/cm 2 , and 
the depletion of ground states and the first excited state 
is only 0.025. Fig. 2(b) exhibits harmonic-generation effi- 
ciencies of at least four orders of magnitude greater than 
that of Fig. 2 (a). It is worth it to be noted that both 
odd and even harmonic peaks appear although the even 
harmonic is weak in Fig. 2, which accords with our above 
analysis that with the increasing internuclear separation 
R, the Rabi splitting of the odd harmonics should con- 
verge towards even harmonics and their intensity should 
become weaker and weaker. Again the calculation result 
also shows the strong-coupling of the ground state and 
the first excited state. 
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Figure 3: Harmonic-generation spectra at uo — 0.05642a. u., 
R = 5.2, I = 6.7 x 10 13 W/cm 2 , initially in the ground state, 
from ID time- dependent exact calculation ((a) and (b)) and 
the formula (14) ((c) and (d)): (a) and (c) the contour of the 
spectrum; (b)and (d) the corresponding fine structure of (a) 
and (c), respectively. The arrows in (a) and (c) show the 
molecular and the atomic like cutoffs of the transfer harmon- 
ics; the arrows in (b)and (d) show the symmetrical splitting 
around the 3th harmonic, and the numbers above the arrows 
show the splitting separation. 



Figure 4: Harmonic-generation spectra at uo — 0. 05642a. it., 
R = 6a.u., I = 5.3 x 10 13 W/cm 2 , initially in the ground 
state, from ID time-dependent exact calculation ((a) and (b)) 
and the formula (14) ((c) and (d)): (a) and (c) the contour of 
the spectrum; (b)and (d) the corresponding fine structure of 
(a) and (c), respectively. The arrows in (a) and (c) show 
the atomiclike cutoffs of the transfer harmonics; the arrows 
in (b)and (d) show the symmetrical splitting around the 3th 
harmonic, and the numbers above the arrows show the split- 
ting separation. 



B. Harmonic spectrum calculated from our model 

Fig. 3 (a) and (c) show the HOHG spectra from the ID 
time-dependent calculation and the numerical calculation 
of the formula (14) with R=5.2a.u. and assuming the 
system initially being in the ground state. Fig. 3(b) and 
Fig. 3(d) show the fine structure of Fig. 1(a) and Fig. 3(c), 
respectively. The splitting separation is Q — XAlw in 
Fig. 3(b) and 1.94m; in Fig. 3(d). The discrepancy is due 
to invalidation of rotating- wave approximation in intense 
field. A plateau in low-order region is also identifiable in 
the spectrum in Fig. 3(c) with a cutoff at the 5th har- 
monic order, which agrees with Fig. 3(a). This corre- 
sponds to the maximum energy acquired by the electron 
in the two- level system in the presence of the field |2?jj. 
Furthermore, a second plateau up to the 25th harmonic 
order which is associated with the ground-continuum 
coupling (referred to as the atomic plateau) can be iden- 
tified. The cutoff position of the atomic plateau can be 
well explained by a semiclassical model|28j. 

Fig. 4 shows the HOHG spectra from ID time- 
dependent numerical calculation and (14) with R=6a.u. 
and E\ - Eq = 0. 0304a. u. The total depletion of the 
ground states and the first excited state is 0.34. The 
splitting separations is Q — 1.2u> in Fig. 4(b) and 1.17a; 



in Fig. 4(d). The position of cutoff in Fig. 4(c) is at 
the 24th harmonic order that agrees with the order 
n = (3.17E7p + I p )/oj, the cutoff law predicted by the 
semiclassical mechanism of HOHG. 

Fig5 shows the comparison of the splitting separa- 
tion(scaled by the frequency u> with u> = 0. 05642a. u.) 
between the model prediction of Q—y/\Q'ji\ + £ 2 (solid 
curve) and the ID time-dependent exact calcu- 
lation (dotted curve) at J2=5.2a.u.(Fig5(a)) and 
i?=6a.u.(Fig5(b)). It is obvious that in the near reso- 
nance case(Fig5(a)), good agreement is obtained while 
the field is weak. However while the field becomes 
stronger, the difference between the theory one and the 
numerical one seems smaller in Fig5(b) than in Fig5(a) 
with the same field intensity. 



C. Influence of initial condition and internuclear 
distance 

At intermediate R, such as R=5.2a.u. or R—6a.u., no 
matter initially the system is in the ground state or the 
first excited state, or the coherent superposition state, no 
distinct difference is found in the spectra from numerical 
calculations. However at large R, the case is rather differ- 
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Figure 5: The splitting separation(scaled by the frequency ui) 
at uj — 0.05642a. u. get from the model prediction of Q(solid 
curve) and from 17? time-dependent exact calculation (dotted 
curve) VS the maximum absolute amplitude of the external 
electric field (atomic unit): (a)R = 5.2; (b) R = 6. 



ent. Fig. 6 shows the comparison of harmonic generation 
spectrum calculated from different initial condition at 
R=8a.u., where E x -E = 0.0052a.u. and Q = 0.95w. In 
Fig. 6(b), where the system initially is in the coherent su- 
perposition state, even harmonics appear although they 
are weaker than the odd harmonics, however in Fig6(a), 
where the system initially is in the ground state, even 
harmonics hardly are distinguishable. The result vali- 
dates our analysis for (20), which predicts that at large 
R if the system initially is in the ground state, only odd 
harmonic can be produced , and if it is in the coherent 
superposition state, even harmonics also should be emit- 
ted although their strength should be weaker than the 
odd ones. If initially the system is in the excited state, 
the case is analogous with the ground state. Other calcu- 
lations at R=16a.u. from different initial conditions also 
show similar result. The position of cutoff in Fig. 6(c) is 
at the 23th harmonic order, which agrees with the order 
n = (3.17f7 p + I p )/u). 

At intermediate R, with the same parameters as in 
Fig4, Fig. 7 shows the comparison of contributions to har- 
monic from the four different terms (10)-(13). It can 
be seen that all of them give equivalent contribution to 
HOHG in this case. 

At large R, with the same parameters as in Fig6(a), 
Fig. 8 shows the comparison of contributions to harmonic 
from the four different terms (16)-(19) if initially the sys- 
tem is in the ground state. It is the term of |0) — ► |p) — > 
|0) that contributes mostly to HOHG in the case. Simi- 
larly, the calculation also shows that if initially the sys- 



Figure 6: Harmonic-generation spectra calculated from ID 
time-dependent exact calculation((a) and (b)) and the for- 
mula (20) ((c) and (d)) at R = Sa.u., I — 6.7 x 10 13 W/cm 2 
and lu = 0. 05642a. u., (a) and (c)initially in the ground state; 
(b) and (d)initially in the coherent superposition state of the 
ground state and the first excited state with equally weighted 
populations. The arrows in (a) and (c) show the atomic-like 
cutoffs of the transfer harmonics. 



tern is in the first excited state, the primary contribution 
to HOHG should come from |1) -> |p) -► |1). 



D. Cutoff law at large R 

Fig. 9 shows the spectra calculated by our analytic the- 
ory and ID numerical simulation with R=16.a.u.. The 
main difference is that there are weak even harmonic 
peaks in Fig. 9(a) but no even harmonic peaks in Fig. 9(a). 
This is because E\ — Eq = has been adopted in the 
analytic formulation but actually these two states are 
not completely degenerated. The position of cutoff in 
Fig. 9(b) is at the 3bth harmonic order, which agrees 
with the order n = (5.33C/ p + ip)/w that is obtained by 
the semiclassical calculation considering complex trajec- 
tories. This is consistent with the finding that at very 
large R, the molecule has different cutoff law 0, El ■ 



IV. CONCLUSIONS 

In summary, with considering the continuum states 
and ionization, the photon-emission spectra of the sym- 
metric diatomic molecule ions in intense laser fields are 
thoroughly investigated and the role of the CR states is 
identified. 
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Figure 7: Harmonic-generation spectra (same parameters as 
in Fig. 4), calculated from the formulas (10)-(13): (a)from 
(10); (b)from (11); (c)from (12); (d)from (13). 



Figure 9: Harmonic-generation spectra at u) — 0. 05642a. u., 
R = 16, / = 1 x 10 14 W/cm 2 , initially in the ground state: (a) 
ID time-dependent exact calculation; (b) calculation accord- 
ing to the formula (20). The arrows in (a) and (c) show the 
atomic-like cutoffs of the transfer harmonics. 
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Figure 8: Harmonic-generation spectra (same parameters as 
in Fig. 6(a)), calculated from the formulas (16)-(19): (a)from 
(16); (b)from (17); (c)from (18); (d)from (19). 



In the near-resonance region {{E\ — E )/uj ~ 1), it 
is shown that the multiple Mollow triplet structure can 
be described by an analytic theory based on rotating- 
wave approximation and the Rabi oscillation between the 
ground state and the first excited state induces the sym- 
metrical splittings of odd harmonics. The splitting sep- 



arati on in weak field can be a pproximately denoted by 
Q=^\ED\ 2 + (E 1 -E -uj) 2 . 

In the strong coupling region ((Ei — Eq)/w ~ 0), all 
of the symmetrical splitting sidebands of odd harmonics 
gradually converge to even harmonics and their ampli- 
tudes get weaker and weaker with the increasing inter- 
nuclear separation R. This is explained by our theory 
in the ultimate condition of {Ex — E )/oj ~ 0(R — ► oo). 
An interesting phenomenon is that in this case if initially 
the system is in the coherent superposition of the ground 
state and the first excited states, there are stronger even 
harmonics emitted than only in the ground or first ex- 
cited state, and their intensity depends on the value of 
\Qr\ = \ED\. 

In all above cases, the molecular plateau and the 
atomic plateau are identifiable in the analytical model 
and the cutoff positions are accordant with the numeri- 
cal simulation and semiclassical calculation. 

We also find that, for the symmetrical diatomic 
molecule, the dipole movement includes four terms, those 
are |0> |p> |0>, |1) |p) - |0>, |0> |p) |1) 
and |1) — > |p) — > |1), respectively, denoting the transi- 
tions from the ground state or the first excited state to 
the continuum state, then back to the ground state or the 
first excited state. It is more complicated than atom case 
because of the existence of the CR states. At intermedi- 
ate R, each process contributes equivalently to molecular 
HOHG. However, at very large R, their contributions de- 
pend on the initial condition. 
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